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Abstract 

We show that when the thermal wavelength is comparable to the spatial 
size of a system, thermodynamic observables like Pressure and Volume have 
quantum fluctuations that cannot be ignored. They are now represented by 
operators; conventional (classical) thermodynamics is no longer applicable. We 
continue the work in earlier papers where quantization rules for thermodynamics 
were developed by analogy with optics and mechanics, by working out explicitly 
the quantum theory of van der Waals gases. We find a wave equation satisfied 
by the thermodynamic wave function as well as solutions (coherent states) that 
are centered at the classical equations of state. The probability of departure 
from the classical theory is dependent on a parameter a which is a property of 
the gas molecule. 
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1 Introduction 



Usually we use quantum mechanics to describe small systems like individual 
atoms and thermodynamics to describe large systems like gases. But there are 
a few examples where quantum mechanics is needed for systems that are large 
enough to be also thermodynamic: a well-established example would be a super- 
conductor. We also now have experiments with quantum gases at temperatures 
small enough that the thermal wavelength is of the order the distance between 
atoms, leading to Bose condensation L lj. Another leap in which the thermal 
wavelengths are comparable to the size of the trap containing the gas can be 
contemplated. That will bring in a new set of quantum effects. 

In this paper we will continue to develop a quantum theory of thermody- 
namics, begun in Ref. [2], By this we do not mean just quantum statistical me- 
chanics. For example, in deriving the partition function of the Bose gas, we use 
quantum mechanics. However, the thermodynamic quantities such as Pressure 
or Volume are then treated as numbers not operators. By quantum thermo- 
dynamics we mea a situation where thermodynamic variables have quantum 
fluctuations and therefore are to be treated as operators. Recall that in classical 
thermodynamics, as in classical mechanics, observables come in conjugate pairs 
such as P, V or T, S. 

In the earlier paper [2] we have suggested a thermodynamic uncertainty 
principle for quantum gases, bounding the product of quantum fluctuations in 
thermodynamically conjugate variables such as P and V. If the gas is cold 
enough, and the trap is small enough, that the thermal wavelength of the parti- 
cles is comparable to the diameter of the trap, we can no longer treat P and V 
as classical (commuting) observables. They are operators acting on a thermo- 
dynamic analogue of the quantum mechanical wavefunction. The equation of 
state of the thermodynamics system is not any more a relation between quanti- 
ties such as P, V,T: such operator relations would be inconsistent. Instead we 
have a differential equation satisfied by a thermodynamic analogue of the wave 
function. 

What is this differential equation, the thermodynamic analogue of the 
Schrodinger equation? The Schrodinger equation of quantum mechanics reduces 
to the Hamilton- Jacobi equation of classical mechanics. Thus a way to discover 
the quantum thermodynamic wave equation is to first find the Hamilton-Jacobi 
formulation of thermodynamics. This was accomplished in Ref. [3]. 

In this paper we will study the resulting quantum wave equation further 
using as an example the van der Waals theory of non-ideal gases. We chose this 
example because it is known to be a good model of real gases: the equations of 
state of hundreds of gases have been fit to van der Waals (vdW) theory, with 
widely available tables [4] of the parameters for each molecule. Moreover, it 
turns out that the vdW is a kind of 'integrable system': by a change of variable 
we discovered in [3] its canonical equations can be solved analytically in closed 

3 Another theory, also called sometimes quantum thermodynamics has already been pro- 
posed in series of interesting papers [5]. Our work does not seem to be related to it . 
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form. Thus this system provides the correct balance between mathematical 
accessibility and physical realism. 

We will work within the canonical rather than the micro-canonical ensemble. 
This means that our approach will not be valid near a phase transition to a 
condensate. A deeper analysis including chemical potential as a thermodynamic 
variable is conceptually similar, but will be more complicated because of the 
extra degree of freedom. We hope to study that extension in a later publication. 

2 The Thermodynamic Uncertainty Principle 

Let us review the argument for an uncertainty principle for thermodynamically 
conjugate variables. Consider a gas of molecules contained by some potential 
well acting as a 'container' or 'trap'. Classically there is a well defined meaning 
to the position of the 'wall' of the container: it is the point at which the particle 
is reflected back into the trap, or the location of the classical turning point. 
But, in quantum mechanics things are not so clear-cut: the particle has a small 
probability to tunnel past the turning point, or to be reflected before the turning 
point is reached. Thus the linear size of the container L has an uncertainty AL 
due to quantum fluctuations, with a resulting uncertainty of order AV = AAL 
in the volume. 

There is also an uncertainty in the Pressure. Pressure is the force exerted 
per unit volume on the wall when particles and reflected back. This is of order 

P = 2ppv. (1) 

Here p is the component of momentum normal to the wall. Also, pv is the 
number of collisions (per unit time per unit area) with the wall: the product of 
number density and average velocity. 

In quantum mechanics, the momentum of the particle and the position have 
quantum fluctuations, with the uncertainties bounded by 

ApAL > | (2) 

This translates to a quantum uncertainty in pressure and L: 

APAL > hpv. (3) 

We should verify that for most conventional applications of thermodynamics, 
this quantum uncertainty can be neglected. To get an order of magnitude of 
this error, we can use the ideal gas equations of state 

PV = Nk B T, (4) 



where ks is Boltzmann's constant and N the total number of particles. Or, 



Thus 

On the other hand, by equipartition of energy, ^kgT is the contribution to 
kinetic energy due to motion normal to the wall; this is equal to the momentum 
times the velocity normal to the wall. So 

where 

Pt = Vmk B T (8) 

We have the equivalent forms 
and since V oc L 3 , 



AP AL V2^A 

P L ~ L U 



where A = -,===&=== is the thermal wavelength of a particle in the gas 171. 

\f2nmKBT or- o i— j 

This is the thermodynamic uncertainty principle. The quantum bose gases in 
magneto-optical traps have achieved very low temperatures already pQ ; it might 
be possible to observe such fluctuations experimentally in a later generation of 
such experiments. 

Thus, if the size of the trap is comparable to the thermal wavelength, 



3V2ttA > L. (11) 

pressure and volume cannot be treated as commuting observables. They are 
quantum operators acting on a thermodynamic analogue of the wave-function 
of quantum mechanics. 

The equation of state therefore cannot be a relation among P, V, T: such an 
operator relation would be inconsistent. Instead we must have a wave equation 
satisfied by the thermodynamic wave function. 

The above condition is to be compared with the much weaker condition at 
which the quantum statistical effects (e.g., Bose statistics) become significant: 

A>(J)* (12) 

That is, the thermal wavelength is greater than the inter-particle distance. 

It is unlikely that a weakly interacting Fermionic system will satisfy the 
condition (fTT|) : the wavelength at the Fermi surface is likely to be too small. 
However, composites of fermions may have this behavior. In any case, weakly 
non-ideal bosonic gases would seem to be the best candidates. 
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3 Equation of State of a Non-Ideal Gas 



For an ideal gas 

where v is volume per particle 



P= T - (13) 



V 

(14) 



k H N 



and ks is the Boltzmann constant^. 

The standard classical model for non-ideal gases is the van der Waals theory. 
It describes very well the behavior of many materials in the real world. Also, 
it gives a good description of the liquid-gas transition away from the critical 
point. 

In the van der Waals theory [6j [7] a molecule is assumed to have a finite 
volume b, so that the actual volume available for the gas is N[v — b]. Moreover, 
there is a force between pairs of molecules. In the bulk the net force is zero as 
there are the same number of molecules in all directions. But at the wall, there 
is a net force normal to the wall (there being no gas molecules on the other side) 
so the pressure is changed by an amount proportional to the square of density. 
If the force between the molecules is attractive this is a decrease in pressure; 
otherwise it is an increase. Thus we are led to a modified equation of state 

^ ' " (16) 



v — b v 2 ' 

the last term being proportional to the density of pairs of particles, or to the 
square of number density. 

Another equation of state is the formula for u = k ^ N , the internal energy 
per particle 

,-fr-i. (.7, 

The last term is the extra energy each particle gains through interactions with 
its neighbors. Thus, we can regard the departure from ideal gas behavior as the 
addition of constants —a, b to the variables v, uv. A thermodynamic description 
in which these are used as the co-ordinates, will be especially convenient. 

In the original van der Waals theory, the molecules attract at large distances 
so that a > 0. But both signs for a are interesting from a modern perspective. 

These arguments apply also to gases at low enough temperatures to be quan- 
tum. For example, in the Lee-Huang- Yang [5] theory of dilute non-ideal Bosc 



4 It will be convenient to scale extensive quantities 

S U 
u ■■ 



? AT k B N 



(15) 



etc. to simplify formulas. Conventional units can be restored at the end using simple dimen- 
sional analysis. 
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gases, the equation of state is of the above form with & = and a < 0: the 
negative sign being because the pairwise force is repulsive. More precisely, 



t2 

a = -8n——a s , (18) 
2mKB 

where a s is the scattering lengtrjf]. If the scattering length is positive the force 
is repulsive. More generally, a is related to the energy of the lowest resonance 
in the two-body elastic scattering of two molecules. We will thus use the van 
der Waals as our generic model of a non-ideal gas even in the quantum regime. 

A gas has a thermodynamic phase space with five co-ordinates (it, T, s, P, v). 
The first and second laws of thermodynamics require the relation 

du = Tds - Pdv (19) 

among the infinitesimal variations of these quantities. Here we see that thermo- 
dynamic variables appear in conjugate pairs such as (T, s) and (P, v). Because 
of the odd-dimensionality, there appears to be an unpaired variable u. However, 
by rewriting the relation as 

du P , 

ds = — + -dv (20) 

for example, we can see that there is nothing special about u as a thermodynamic 
co-ordinate. As explained elsewhere [31 HO] classical thermodynamics has a 
natural formulation in terms of contact geometry, the odd dimensional analogue 
of the symplectic geometry that describes mechanics. The quantities (u, v) are 
canonically conjugate to 

IP 

Pa = ji, p v =j t - (21) 

Among the five variables there are three relations, the equations of state of 
the gas. These relations are most simply described by giving the entropy s as a 
function of (u,v): the fundamental relation. For the van der Waals gases, this 
function also depends on two additional parameters (a, b) which describe the 
departure of the gas from being ideal: 

s = s(u, v; a, b) (22) 

Given this function, the remaining variables are given by the derivatives 



T ~ (dy,).' T~(dv 



(23) 



For fixed (a, b), this determines a two dimensional surface in the five-dimensional 
thermodynamic phase space: the Lagrange sub-manifold of this particular gas. 



5 fes appears because our definition of energy includes it as a factor: u = jrr—- We have 
to convert from the units used in Ref. |8] to ours. 
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The fundamental relation of case of a van der Waals gas is, explicitly, 

3 3 

s = log(v — b) — — log v + — loghw + a] . (24) 

A simple calculation of the derivatives shows that we get the van der Waals 
formula for P upon eliminating u. 



4 Hamilton- Jacobi Theory of van der Waals Gases 

It is possible to 'unify' [3] the van der Waals model of gases, so that the funda- 
mental relation s(u, v; a, 6) is the solution of a common differential equation: the 
parameters a, b appear as constants of integration. This is the analogue, in ther- 
modynamics, of the Hamilton-Jacobi formulation of mechanics or optics. In the 
usual picture where u, v are used as co-ordinates, entropy would the analogue 
of the eikonal. 

By eliminating a, b we can get a relation among the five thermodynamical 
co-ordinates: 

F(s,u,v,p u ,p v ) = (25) 
If we make the canonical substitutions 

ds ds 
Pu= ^? Pv= d~v (26) 
we get the thermodynamic analogue of the Hamilton-Jacobi equation: 

ds ds, 
du' dv' 

All solutions of such a first order PDE are determined by a fundamental 
solutiorH that depends on two constants of integration. It is in some ways 
analogous to the plane wave solution of a wave equation; every other solutions 
can be expressed as a linear superposition of plane waves. All other solutions of 
a first order PDE are obtained from the fundamental solution by the steepest 
descent approximation to such a superposition: by an extremization. 

The explicit form of the equation for the hypersurface of the vdW famility 
of gases (obtained by elininating a, b from the equations of state) is [3]: 



F(s,u,v, 



) = o 



(27) 



n 2 



Vp v - Up u + 



pi 
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v 2 e 



so that the H-J equation is 

'd£ 
dv 



ds 
du 



ds_ 

du 



27 



v 2 e~ 



2 s 



(28) 



(29) 



We will now solve the Hamilton-Jacobi equation and recover the fundamental 
relation (f2"4"]l . This exercise will also help us identify the natural variables of the 
system, in which the quantization will be easier. 



6 See where it is called the 'complete integral'. We propose the name 'fundamental 
solution' instead, which fits better with the term 'fundamental relation' as used in thermody- 
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5 Normal Co-ordinates 



It will be convenient to change to some new variables (rather like normal co- 
ordinates of a mechanical system) (</>, qi,q 2 ,Pi,p 2 ) which still satisfy thermody- 
namic relation 

d(j) = pidqi + p 2 dq 2 . (30) 

That is, a Legendre Transformation, the thermodynamic analogue of a canonical 
transformation of classical mechanics. We can choose the transformation such 
that the H-J equation becomes independent of 4>- This will make it possible 
to find the fundamental solution by separation of variables and recover the van 
der Waals relation. These variables were discovered in [3], so we just give the 
answer and show how the reader can verify that it has the correct properties: 



qi = «, 92 = ^ uv > ( 31 ) 



and 



Pi 



3 

Vp v - Up u + - 



[v e 



P2 = \pu{v~ 1 e s Y , (32) 



<P=^vhi s (33) 

It is straightforward (and not too tedious) to check that this is indeed a Legendre 
Transformation: 

d(f) = pidqi + p 2 dq 2 ds = p u du + p v dv. (34) 



These variables have been chosen such that, the hypersurface equation (125|) 
is simply 

p\pI = 1- (35) 
The H-J equation is then obtained by the canonical replacement 

d<P d<j> 
Pi = «— , P2 = -k— (36) 
dqi dq 2 

to get 

We can now solve this PDE by separation of variables, and recover the 
equations of state of the vdW gas: 

<f>(qi,Q2) = <Pi(qi)<t>(<&) (38) 

so that 

[0'i(9i)] 2 [02(g 2 )] 2 [0 2 fe)] 3 [0i(g 2 )] 3 - 1. (39) 

Or, 

$ 2 =C 2 <f 2 l (40) 



8 



for some separation constants satisfying 

C?Cf = 1. 
These ODE are easy enough to solve: 

</>i(<?i) = Ki(qi ~ ^2(92) = ^2(92 - q' 2 Y l 



(41) 



(42) 



Here, q[,q' 2 are constants of integration. The product K\K 2 is determined 
through 

2 s 3 5- 

Ci = -K?, C 2 = -Ki. 



so that 

Thus 

where 



5 x ' ' 5 
2 2 3 3 [K 1 K 2 ] 5 = 5 5 =>■ ^1X2 = 5 2"? 3~s 

0(91,92) = - <7i)^(<?2 - 9 2 )^ 



Recalling the definitions d> = |w = e5 s ,gi = t>, ga — l uu we get the familiar 



(2 2 3 3 ) 

formula (|24[) for the entropjjj of a vdW gas: 

s = log 



(43) 

(44) 
(45) 
(46) 



3"| 

(v-b)[u + - 

v 



with the identification of the constants of integration: 



b = q[, a = --q 2 . 



(47) 



(48) 



Thus our Hamilton-Jacobi equation 'unifies' the vdW gases: the fundamental 
relation of different gases are solutions of the same PDE, with the parameters 
describing the gases appearing as constants of integration. 
Using pi = g^,P2 = J^j) we get the equations of state: 



Pi 



2 q 2 - q' 2 



3 qi-q'i 



P2 



2 q 2 -q' 2 



3 qi - q[ 



(49) 



6 The Wave Equation for a Quantum vdW Gas 



We can now write down a Schrodinger equation by the rule 

d 

Pk -> -i^r~ 



(50) 



7 Up to an additive numerical constant. 
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Pi 



6 
5| 
4| 




Figure 1: The curve p\p\ = 1- The wavefunctions we study are supported on 
the branch with p\ > 0. 



familiar from quantum mechanics and optical 

d 5 



dqfdq? 



V> = V> (51) 



Recall that the equation we are solving is invariant under translations in 
9i, <?2- Using this symmetry, we can get the Fourier representation 

^( q ) = [ e ^i+P^} S{p 2 p 3 _ 1)x{pup2)dpidp2 (52) 



with an arbitrary function x on the curve (analogous to the 'mass- shell' of the 
Klein- Gordon equation) 

p\pl = !• (53) 

This curve always has p2 > and has two branches with p\ > and p\ < 0. 
If we look back at the definition of p\ we will see that positive pressure implies 
Pi > 0, which is the branch we will choose. Note also the scaling symmetry 

Pi -> A 3 pi, p 2 -> A" 2 p 2 , A > 0. (54) 

This is analogous to Lorentz transformations that leave the wave equation un- 
changed. A parametric solution in terms of rational functions of the curve is 

Pi=k 3 , P2 = k- 2 . (55) 

Any solution of the wave equation supported on the positive branch of the 
shell above is of the form 



8 We could easily transform to the more familiar thermodynamic variables through qi = 
v > 12 = 3 uv but then, the equation will look quite complicated. Also, it does not help in 
solving the equation. So we will continue to work with the 'normal co-ordinates' q\,q2- 
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= ^ e^ k3+ "^ 2 }m^. (56) 



We choose ^ as the measure of integration to respect the scaling symmetry 
k -* Xk. 



7 States and Observables 

The quantum Hilbert of the system is the space of solutions to the wave equation. 
Thus it can be identified with the space of functions ip(k) on the shell. It seems 
reasonable to postulate that the inner product should be 

\m\ 2 ^. (57) 

K 

We have chosen the norm such that the scaling 

Mk) = ^P(Xk) (58) 

leaves the norm invariant : 

ll^|| 2 = |M| 2 . (59) 

In this picture, the 'pressure' variables (analogous to 'momentum' variables 
of quantum mechanics) are just multiplication operators. 

Pi=k 3 , P2 = k- 2 . (60) 

But the 'co-ordinate' variables <7i,<?2 will be differential operators [12] . Up to 
additive constants, q± = igp-,52 = ~i gjtr=i ■ Taking into account the condition 
that they be hermitearj^l w.r.t.the above inner product, we havJ^l 



qi=q i + 6 



k 3,k §k 



i 1 
<l2 = q 2 + - 



k 2 k — 



(61) 



The quantities ^,Q^ are constants that don't affect the commutators of the 
operators. 

The expectation value of an observable A represented as a differential oper- 
ator in k is then 

f?\mn " m 

We now look for a state that is of finite norm and is centered at the clas- 
sical solution. More precisely, the expectation values of the operators above 
Pi>P2, Qi, 92 must satisfy the classical equations of state. 

Note that in this interpretation, the co-ordinate space wavefunction ip{qi, 92) 
does not have the meaning of a probability amplitude: in fact it is usually not 



9 Noice that in this picture, i-J^ is not hermitean, but ik-^ 
10 [A, B]+ = AB + BA is the anti-commutator 
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Figure 2: The wave function with x = 0, k — 5,(7=1 



even square integrable. This is similar to the case with Klein-Gordon theory 
PI- 

There could be many such states, but the simplest would be 'plane wave' 
solutions, which are the wave analogues of the fundamental solution of the 
Hamilton- Jacobi equation. Any other solution of the wave equation would be 
super-positions of plane waves. For these 

$(fc) = 0(k)e l[zik3+Z2k ~ 2] (63) 

for constants!^ zi . z%. But, as in quantum mechanics, these are not normahzable 
if zi,Z2 are real. If they have a positive imaginary part, they become square- 
integrable: 

Im zi,Im z 2 > 0. (64) 

The parameters z\ , z 2 can be chosen so that the mean values of observables are 
centered at classical values: loosely analogous to the coherent states of quantum 
mechanics [13] . 

We are thus led to the wave function 

4>(k\K, x, a) = e ~ lxk3 e~^ & iT+H^Y] ( 65 ) 



This function is smooth: it vanishes at k — along with all of its derivatives. 
The minimum of the quantity in square brackets is at k: this is the most probable 

11 The Heaviside function in the expression of tp is actually not necessary, because by defi- 
nition of the state space, the wave function is defined only for positive values of k. 
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0.002 - 
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-0.002 - 
-0.004 - 
-0.006 - 



Figure 3: The real and imaginary parts of the wave function for x = 1,k 
5,a = l 

value of k. Near that point, a little calculus gives, 

(fc-«) 2 



e 2 ° 



(66) 



In the limit of small a, we can estimate expectation values using the steepest 
descent approximation; essentially, approximate it by the Gaussian: 



< k k, Ak 2 =< [k— < k >) 2 a 2 



It follows that 
Also, 



< pi k 3 , < p 2 >« k 2 . 



Apx w 3k 2 ct, Ap 2 ~ 2kT 3 (7. 
In the same approximation, 

3 

< <?i >« q[ + X, <q 2 >« q' 2 + ^ 5 X. 

Eliminating x, k from the formula for the expectation values, 

_ T2 < g 2 > -ga l * 
[3 <qi> -q[_ 

we get the classical equations of state (compare with equation |49| 
~2<q 2 >-q' 2 ]% \2<q 2 >-q' 2 



< Pi >t 



i<qi> -q[ 



3 < qi > -q[ 



(67) 

(68) 
(69) 

(70) 
(71) 

(72) 
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A little more work shows that 



6k 2 



1 




(73) 



Thus, for small <r, 



ApiAqi w -, Ap 2 Aq 2 w - 



(74) 



Although the probability of small fluctuation is well-approximated by a Gaus- 
sian, rare large fluctuations can have a very different probability. If the strategy 
for observing such fluctuations rely on such rare events, the full wave-function 
given above must be used and not just its Gaussian approximation. 

Thus we have found a state that is centered at the classical thermodynamic 
values of the observables, with fluctuations determined by the parameter a. We 
expect that this parameter is a property of the ground state of the molecule 
that makes up the gas, just like the van der Walls constants a, b. For example b 
is the volume of the molecule: b~s is roughly the mean distance of the outermost 
electron from the center of mass of the molecule. Similarly |, which has the 
dimensions of energy is roughly the position the lowest lying resonance in two 
body elastic scattering of the molecules. Just as van der Waals theory itself 
makes no prediction for a, 6, we make no prediction within our theory for a. 
However, we should expect that it is roughly the same order of magnitude as 
a, b. For the reasons explained earlier, quantum fluctuations in thermodynamic 
quantities are suppressed when the thermal wave length is small, so even when 
a is of the same order as a, b, its effect will be measurable only at very small 
temperatures. 

Although a, b used to be thought of as intrinsic properties of molecules, it is 
possible to manipulate the wavefunction of the outermost electron by external 
electromagnetic fields. Feshbach P3] resonances can be used to change the 
scattering lengths by several orders of magnitude, leading to the same change 
in a. Perhaps such clever experimental techniques will be found some day to 
enhance the value of a as well. 

It would be most interesting to measure experimentally the quantum fluc- 
tuations in thermodynamic quantities. In the theoretical direction, the analogy 
with Klein-Gordon theory suggests a second-quantization of our wave equation. 
We do not know yet what that would mean physically. 
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9 Appendix: The van der Waals Constants 

Tables of van der Waals constants [4] are easily available. We recall a few values 
to give an idea of the order of magnitudes. 
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For hydrogen molecules, b — 0.027 x 10 3 ^j- Using Avogadro's number, 
N A = 6 x 10 23 we get 

6= 0^ xl0 -26 m 3 M ( 3 _ 56xl0 -lO m )3 (?5) 

Thus b is of the order of the volume of a cube of side equal to 7.1 Bohr radii. 

Also, for hydrogen molecules, a = 0.025^p. Thus | w 1^ w lO^eV 
which is much smaller than atomic energy levels. 

For Helium, b = 0.024 x lO" 3 ^, so that &i w 6.8a B and a = 0.0034^, so 
that | « 1.5/xeV. As expected, Helium atoms interact much more weakly than 
Hydrogen molecules. 
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